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Abstract 

A model of diluted random field sustained by quenched volume de- 
formations is shown to reproduce puzzling physical features found in 
X(CN) x Yi_ x mixed cyanide crystals. X is an alkali metal (K, Na or 
Rb) and Y is a spherical halogen ion (Br, CI or /). Critical thresholds 
x c at which associated first order ferroelastic transitions disappear are 
calculated exactly. The diluted random field is shown to compete with 
compressibility in making the transition first order. Transitions are then 
found to remain first order down to x c except in the case of bromine 
dilution where they become continuous. All the results are in excellent 
agreement with available experimental data. 
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1 Introduction 



In recent years a great deal of experimental effort has been devoted to inves- 
tigating the existence of a possible orientational glass state in mixed cyanide 
crystals [1, 2]. In particular diluted cyanide mixtures X(CN) X Y]_- X where X 
is an alkali metal (K, Na or Kb) and Y stands for a spherical halogen ion (Br, 
CI or I) are known to exhibit a series of puzzling physical features. 

For instance, the ferroelastic transition which occurs in pure alkali-cyanides 
XCN [3] is found to disappear at a cyanide concentration threshold x c which 
is not equal to the underlying lattice percolation threshold. The actual value 
of x c varies with X and Y compounds. In the case X = K it is x c = 0.60 and 
x c = 0.80 respectively for Y = Br and Y = CI. 

In addition, below x c there exists experimental evidence for local freezing 
of cyanide orientations, together with the absence of quadrupolar long range 
order. In parallel X(CN) X Y\- X mixtures remain cubic at low temperatures [1]. 
Analogies to spin-glasses and random field systems have thus generated much 
interest, in particular in elucidating the physical nature of the orientational 
glass state [1, 2]. 

Moreover the physics above x c is still puzzling. The transition is found 
to be continuous at x c for X(C N) x Br\_ x [4] while it stays first order with 
X(CN) x Ch- x [5]. 

In this paper a model of diluted random field sustained by quenched volume 
deformations is proposed. It is shown to reproduce most of puzzling physical 
features found in X(CN) x Yi- x mixed cyanide crystals. 

A derivation of a microscopic calculation of critical thresholds x c is pre- 
sented. On this basis, an Hamiltonian of a diluted compressible ferromagnetic 
Ising system is studied within mean field theory. Diluted staggered fields are 
then included to account for the quenched volume deformations generated by 
the substitution of spherical halogens to dumbbell-shaped cyanides. 

Compressibility is found to produce a first order transition at cyanide con- 
centrations larger that some threshold Xl- By contrast, diluted random fields 
are shown to activate a first order transition only below another cyanide con- 
centration x r . 

Transition orders are then discussed for a large variety of mixtures. Rela- 
tive amplitudes in x c , xl and x r are instrumental in determining the transition 
order. Transitions are predicted to remain first order down to x c for all mix- 
tures except in the case of dilution with bromine which exhibits a continuous 
transition at x c . Our results fit perfectly to available experimental data. 

2 More than site dilution 

Starting from pure XCN, substitution of one spherical molecule Y to a dumbbell- 
shape cyanide ion has first the effect of diluting one rotational degree of free- 
dom. Within a simple site dilution picture, associated phase diagram would 
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be only slightly modified by dilution. The transition temperature would be a 
decreasing function of the "hole" (no degree of freedom) density. Moreover it 
would vanish at the lattice site percolation threshold. Here, it would mean fer- 
roelastic transitions disappear at x c = p c for all X(CN) X Y]_- X mixtures where 
cyanides are on a fee lattice where p c = 0.198 [6]. However, from experiments 
%c 7^ Pc and x c varies from one system to another. 

Several physical arguments can explain the discrepancy between x c and 
p c . Most of them like for instance the existence of long range interactions 
would lead one to give up the site dilution approach. However it could be kept 
allowing simultaneously an effective density of free-to-reorient cyanides. In 
other words, dilution would have the additional effect of "neutralizing" some 
cyanides taking them off the process of propagation of quadrupolar long range 
order active in the ferroelastic transition at stake. 

The above picture is given a physical basis by postulating the existence of 
an effective density of free-to-reorient cyanides Xf, such that all ferroelastic 
transitions disappear at the same threshold Xf jC = p c . It implies in parallel a 
density of "neutralized" cyanides x n which depends on the mixture itself with 

At this stage we do not need to elaborate on the physics of these "neutral- 
ized" cyanides. What matters here is their non-participation in the propaga- 
tion of long range order. The effect being a function of the mixture, we define 
naturally 

x n = a(l-x), (1) 
where a is the average number of cyanide "neutralized" by one Y ion. 



3 A geometrical model 

To implement our description, we need a quantitative calculation of a. Notic- 
ing XCN in bulk is cubic only on average, a local deformation of the unit 
cell is expected while substituting one spherical Y ion to the dumbbell-shaped 
cyanide. Let us then assume local XY unit cell embedded in XCN, to match 
that of XY in bulk. Such a mechanism yields a local volume deformation of 
former unit cell of amplitude, 

&v = ^(a 3 KCN -a 3 KX ) , (2) 

(4 molecules per unit cell). Each one of the c F-nearest neighbours (c = 12 on 
the fee cyanide sublattice) being affected by these local volume deformations, 
the overall volume deformation affects on average a total volume AV = cAv. 

Volume deformation AV in turn induces some deformation of neighbouring 
unit cells. The deformable part of these cells is free volume which originates 
from volume differences between a given molecule in bulk vxcn and actual 
molecule itself v a . We thus have a free volume per molecule, 

Vf = vxcn - v a = — {v x + v C n) , (3) 
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where clxcn is the pure XCN lattice constant, vx and vqn are ion volumes 
[7]- 

Dumbell-shaped CN reorientations are directly coupled to the unit cell 
shape via steric hindrance mechanisms. The existence of a finite number of 
equivalent orientations results from the difference in symmetries between re- 
spectively the molecule and the lattice. Symmetry being lower for the molecule, 
any volume deformation of a cyanide cage (increase or decrease) lowers the cor- 
responding symmetry which in turn disturbs the regular dynamics of cyanide 
reorientations. More precisely the number of accessible orientations is lowered. 
Here we are assuming that a cyanide whose cage is deformed (increased or de- 
creased) does not participate in the propagation of long range order at stake 
in the ferroelastic transition. Accordingly one Y substitution will affect on 
average, 

a = c- , (4) 

v f 

unit cells by deforming their respective free volumes. Along these lines, the 
effective cyanide density of free-to-reorient cyanides is, 

Xf = x — a{l — x) . (5) 

The equality Xf >c = p c yields the critical threshold x c , 

* = ^. (6) 
1 + a 

Below x c a region of randomly oriented ferroelastic domains with no static 
phase transition is obtained. These domains shrink with increasing dilution to 
disappear eventually when Xf = at a new threshold [8], 

1 + a 

Comparing Eqs. (6, 7) shows it is the non-zero value of p c which produces 
the randomly oriented ferroelastic domains. Domains are thus created by site 
statistical fluctuations. Shear torque experiments [9], as well as diffraction 
experiments [10] suggested a very similar phase diagram. 

Without a fitting parameter direct microscopic calculation of x c and Xd is 
readily performed using cristallographic data [11, 12]. The cyanide ion volume 
is vcn = 34.47A 3 [13]. The results are obtained for various mixtures (see the 
Table) including systems for which no experimental data are available. 



4 Dilute random fields 

Within the simplest model, cyanide orientations can be represented by ferro- 
magnetic Ising spins {Si}. A random site variable q is then introduced to 
account for site dilution. It is one if site % is occupied by a cyanide and zero 
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otherwise. We have {ei} av = x, where {...} a v denotes a configurational average 
over site disorder. 

Elastic degrees of freedom are also introduced via volume fluctuations to 
account for the first order character of pure XCN transition. Within an 
harmonic model they are integretated out to yield an effective Hamiltonian, 

H eff = ~G e i e jSiSj - E i XI e i e jSiSj) 2 , ( 8 ) 

<i,j> <i,j> 

where G and E are constants [14]. 

With respect to the "neutralization" phenomenon two mechanisms seem 
possible to block some dipoles in a Y ion vicinity. It could be either a local 
increase in orientational barriers or the existence of quenched orientational 
random fields. Local symmetry is the major difference between them. Cubic 
symmetry is preserved in the first case and not in last. However global cubic 
symmetry is preserved for both. Available experimental data do not allow 
actual freezing process to be selected. 

Within a two equivalent orientation model, a cyanide prevented from reori- 
entation is indeed trapped along one direction. On this basis, probability p t to 
have a local random field is identical to the probability of having a deformed 
cyanide cage. From the density x t of trapped cyanide we obtain, 

Pt = S a f^ if*-<s<l . (9) 

[1 if X < Xd 

Associated contribution to the Hamiltonian is, 

# t = -E e ^, (io) 

i 

where hi is a quenched diluted non-symmetry breaking random field. The 
probability distribution function is, 

P(fn) = ^[5(h t -h) + 5(h t + h)] + (1 - Pt )5(h t ) , (11) 

which satisfies both required symmetry conditions P{hj) = P(—hi) and hi = 
where the overline denotes configurational average over the fields. Adding this 
term to H e ff results in full Hamiltonian H T = P e jj + H t which is rather 
complicated. 

A partial mean field treatment is now performed (more details will be 
published elsewhere) to extract main physical features of the model. The 
order parameter is m = {< tiSi >} av where averages are done over both site 
and random field disorders. The double configurational averaged lattice site 
free energy is, 

T = -cGm 2 H m 4 

2 4 

Pt 



-{ln[cosh(f3cGm + f3Lm 3 + f3h)\ + ln[cosh((3cGm + (3Lm 3 - (3h)}} 
2 



— xIcbT 

+ (1 - p t )ln[cosh(l3cGm + (3Lm s )} - k B Tln(2) , (12) 
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where (3 = -r^, &b is the Boltzman constant and T is the temperature. Site 
and RF disorders have been assumed to be uncorrelated in the averaging, i.e., 
variables q and hi are statistically independent. The actual value hi = h or 
—h is independent of the value of e« at a given site i. Also, whether hi = ±h or 
does not depend on site i itself. A statistical dependence between hi and e« 
could only results from spatial correlations between CN sites which are anyhow 
neglected in a MFT. 



5 Tricritical points 

An exhaustive study of Eq. (12) is out the scope of the present work. At 
this stage it is useful to analyse two simple limit cases which are physically 
meaningful. First it is the zero-steric hindrance effect case with a — 0. It 
implies x c = p c , Xd = and p t — 0. 

From a Landau expansion of Eq. (12) a continuous transition occurs at the 
critical temperature, ksT c = xcG under the condition of a positive quadratic 
coefficient, 

B = ~h + h< < 13 > 

which results in the condition x < xl on cyanide concentration where, 

xl = C§) 1/2 ■ (14) 

At x = xl the continuous transition along k B T c = xcG turns first order via a 
tricritical point (B = with a positive free energy sixth order coefficient). 

At x — 1, pt — (Eq. (9)) even if a ^ 0. From experiment, pure XCN 
exhibits a first order transition. Therefore all plastic systems must satisfy 
x L < 1 which sets condition L > Dilution is then found to weaken the 
first order character of the transition. Associated negative quartic term B gets 
smaller in amplitude to vanish eventually at a tricritical point (x — x£)- There 
the transition is continuous with tricritical exponents. Upon further dilution 
the transition becomes second-order for < x < xl. 



6 Results 

We now consider the zero-compressibility case (L = 0) to study random field 
effects (a ^ 0). At x — 1 (no dilution) with p t an independent external 
parameter, the zero-compressibility free energy becomes identical to that of 
trimodal random field Ising models [15, 16, 17] which have been previously 
studied. A first order transition is found only for 0.73 < p t < 1 and for some 
restricted range of random field intensities ~ 0.55 < ^ <~ 0.65 [17]. 

These results produce an additional threshold in cyanide density around 
p t = 0.73 (Eq. (9)), 

Xr " otsT^ • (15) 
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Only at x < x r dilution can turn the transition first order via random fields. 
However condition x c < x r must also be satisfied since the transition itself 
disappears at x c . The equivalent constraint on a gives, 

Eq. (22) gives x r = 0.60 and x r = 0.83 for respectively K(CN) x Bri_ x 
(a = 1.12) and K(CN) x Ch- x (a = 3.48). In parallel we found x c = 0.62 
and x c = 0.82 for Br and CI dilution respectively (see the Table). On this 
basis we conclude that upon dilution, while random fields can turn the tran- 
sition first order in K{CN) x Cl\- x mixtures (x c < x r ), they cannot do it for 
K(CN) x Br\- x mixtures (x c > x r ). 



7 Discussion 

It is worth noting that in the vicinity of tricritical points mean field results 
are indeed reasonable. However, at this stage it is of importance to stress that 
some doubt exists with respect to the validity of the mean field result of a 
first order transition for the random field Ising model. A bimodal distribution 
yields a tricritical point while a Gaussian distribution does not [15, 18, 19, 20]. 

In parallel Monte Carlo simulations showed a tricritical point to occur in 
diluted ferromagnets in a staggered field at weak dilution once next nearest 
interactions are introduced [21]. An exact mapping was earlier found out 
between the mean field treatment of respectively a bimodal random field and 
a staggered field [22]. 

Since our results are in agreement with experiments we have to validate our 
mean field treatment of a bimodal distribution. This is achieved introducing a 
staggered symmetry among diluted random fields. It seems to indicate volume 
deformations may produce staggered random fields. 



8 Conclusion 

In summary, compressibility and staggered random field are active in making 
the transition first order. However compressibility weakens with dilution to 
fade out in the vicinity of x c for all systems which means it is reasonable 
to assume Xl < x c . On the opposite staggered random fields start to be 
instrumental at dilution below x r which has an effect only when x r > x c , i.e., 
a > 2.09 (Eq. (16)). From the Table this condition is found to be always 
satisfied except for mixtures with bromine. 

We can thus predict dilution with chlorine and iodine keeps the transition 
first order down to x c . Only bromine turns the transition continuous. A 
tricritical point is thus expected for X{CN) x Br\_ x mixtures. At this stage our 
predictions reproduce experimental results with respect to potassium systems. 
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Additional experiments on mixtures with sodium and rubidium will check the 
validity of our model. 
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Table 1: Numerical values calculated for a, x c and Xd and experimental thresholds 
when known (denoted by 'exp:'). Units for lengths and volumes are A and A 3 . 
See details in the text. Error bars are within respectively ±0.01 for all data in the 
Table, and ±0.05 for experimental thresholds. 



XCN/Y 


Q-xcn 


axY 


Av 


v f 


a 




Xd 


KCN/Cl 


6.53 


6.29 


7.30 


25.19 


3.48 


0.82 exp: 


0.80 


0.78 exp: 0.75 


KCN/Br 


6.57 


6.60 


2.36 


25.19 


1.12 


0.62 exp: 


0.60 


0.53 exp: 0.50 


KCN/KI 


6.53 


7.06 


18.46 


25.19 


8.79 


0.92 exp: 


0.90 


0.90 


NaCN/Cl 


5.90 


5.65 


5.99 


12.67 


5.68 


0.88 exp: 


~ 0.80 


0.86 


NaCN/Br 


5.90 


5.97 


2.11 


12.67 


1.99 


0.73 




0.67 


NaCN/KI 


5.90 


6.47 


16.63 


12.67 


15.74 


0.95 




0.94 


RbCN/Cl 


6.82 


6.58 


8.05 


30.98 


3.12 


0.80 




0.76 


RbCN/Br 


6.82 


6.85 


1.19 


30.98 


0.46 


0.45 exp: 


0.55 


0.32 
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